Prom perfect to fractal transmission in spin chains 
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Perfect state transfer is possible in modulated spin chains [Phys. Rev. Lett. 92, 187902 (2004)], 
imperfections however are likely to corrupt the state transfer. We study the robustness of this 
quantum communication protocol in the presence of disorder both in the exchange couplings between 
the spins and in the local magnetic field. The degradation of the fidelity can be suitably expressed, 
as a function of the level of imperfection and the length of the chain, in a scaling form. In addition 
the time signal of fidelity becomes fractal. We further characterize the state transfer by analyzing 
the spectral properties of the Hamiltonian of the spin chain. 
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I. INTRODUCTION 

The ability to transfer a quantum state between dis- 
tant parties is one of the basic requirements in many 
quantum information protocols, we mention for example 
quantum key distribution or teleportation 0] . A very 
successful area where the implementation of quantum 
state transmission has been realized is quantum optics. 
The carriers of information (photons) can be addressed 
and transmitted with high control and with a low level 
of decoherence. Very recently, in view of the great po- 
tentialities of solid-state quantum information, attention 
is also focusing on the problem of the transfer of quan- 
tum information in a solid-state environment. A possible 
way to follow would be to properly design couplings be- 
tween optical and solid-state systems Alternatively 
one could also think to realize quantum channels using 
condensed-matter systems. In Ref. Q Bose has shown 
that a Heisenberg spin chain is able to act as a quantum 
channel over reasonable distance (~ 10 2 lattice sites). 
Information capacities for this Heisenberg channel have 
been analyzed in Ref. In Ref. [|| a slightly differ- 
ent scheme has been proposed, in which the simple spin 
chain has been replaced with an isotropic antiferromag- 
nctic spin ladder. A great advantage of these approaches 
is that state transfer occurs due to the interaction be- 
tween the spins of the chain and no dynamical control is 
required (except for the preparation and the detection of 
the state). Proposals to implement this scheme with su- 
perconducting nanocircuits 0, @ have been already put 
forward and very likely these implementations can also 
be extended to other solid-state systems. 

Perfect transfer can be achieved over arbitrary dis- 
tances in spin chains under many different hypotheses: 
by a proper choice of the modulation of the coupling 
strengths as suggested || , if local measurements on the 
individual spins can be implemented [Icj . when commu- 
nicating parties have access to limited numbers of qubits 
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m a spin ring [11| or by using several spin chains in 
parallel |12|. As in all situations in quantum informa- 
tion, the efficiency of such protocols relies on the capa- 
bility of isolating the experimental setup from the ex- 
ternal world (decoherence) , and on the possibility to re- 
duce all possible static imperfections [13, Jl4l ] . In the 
case of the protocols presented in Refs. , where no 
control is needed on the system during the state trans- 
fer, that is, no dynamical control is applied, the cou- 
pling to the environment is supposed to be weak. It 
thus remains to be seen how the quantum channel is 
robust against static imperfections which would be un- 
avoidable (especially in solid-state implementations with 
engineered nanodevices) . Particularly important to con- 
sider is the case when, in the ideal case, perfect state 
transfer is obtained. 

In this paper we address the problem of the effects of 
static imperfections on the protocol presented in Ref. , 
which has already been experimentally implemented in 
Ref. |15| using a three qubit nuclear magnetic-resonance 
quantum computer. We study the sensitivity of the state 
transfer to random variations both of the coupling be- 
tween the spins and of an externally applied magnetic 
field. In view of the possible applications with solid- 
state systems, also the case of correlated disorder will be 
considered. Similar questions for quantum computation 
protocols have been already analyzed in Refs. |l3l Il4| . 
In that case the loss of efficiency of the protocol was 
related to the appearance of quantum chaos in a quan- 
tum computer register. This relation has been charac- 
terized studying the level spacing statistics. Following 
these lines we study the transition of the level spacing 
statistics of the spin chain in the presence of static im- 
perfections. Even though it is not possible to frame this 
problem with the random matrix theory [16|, we show 
that the level spacing statistic is still a convenient tool 
to describe the system efficiency in performing the state 
transfer. 

The presence of static imperfections leads to another 
clear signature of the modified properties of the spectrum 
in the fidelity. The degradation of the state transfer cor- 
responds to the emergence of a fractal signal, i.e., the 
fidelity changes from a periodic function of time to a 
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fractal time series. This behavior has the same origin as 
the one found in the probability densities of the quantum 
evolution in tight-binding lattices 0, 0] . 

The paper is organized as follows. In Sec. [n] we in- 
troduce the model used throughout this work, we set up 
the notations, and we briefly review the quantum state 
transmission protocol of Refs. 0, ||. We then analyze 
the fidelity of the transferred state (Sec. Hlljl . and the 
level spacing statistics of the Hamiltonian of the system 
(Sec. HVfl . Interestingly, the dependence of the fidelity, 
as a function of the length of the chain and the level of 
disorder, obeys simple scaling laws. In Sec. we take a 
closer look at the behavior of the fidelity as a function 
of time. The presence of static imperfections leads to 
a fractal behavior of the time signal of the fidelity. In 
the same section, we relate the fractal dimension to the 
amount of disorder present in the chain. The last section 
is devoted to the conclusions. 



II. MODEL 

The protocol introduced in Ref.Q enables quantum 
state transfer between two parties by means of a spin 
chain: The state of the left-most qubit is transferred to 
the right-most qubit after a given time (dictated by the 
dynamics of the chain). In Ref.||J the approach is the 
same as in Ref.0, the idea is to use a modulated chain 
whose Hamiltonian is given by 
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In Eq. (1} iV is the number of spins in the chain, a 
<t|, cr| are the Pauli operators of the kth spin. The pa- 
rameters Bk and J k are, respectively, the local magnetic 
field and the exchange coupling constant. Both couplings 
depend on the position of the site (or of the link) on the 
chain. In order to achieve perfect state transfer, the sys- 
tem parameters are chosen to be 



B k =0, J k = Jy/k(N - k) . 

The spin chain is initially (at time t = 0) prepared in the 
state 



|#o(tf,¥>)> = (cost?|0)+sin^e^|l))«)|0) c 
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that is, the left-most spin is prepared in a given superpo- 
sition of its two levels while the others are in their ground 
state. The state J2Jl will evolve accordingly to the dynam- 
ics dictated by Eq. ([Q). Since the Hamiltonian commutes 
with the total spin component along the z direction, the 
relevant sector of the Hilbert space is spanned by the 
states 




FIG. 1: Fidelity of the Nth spin with N = 100 as a function 
of time, a) Without imperfections, e.j — eb — 0. b) With 
imperfections ej = 10 -2 , eb = 0, N a v = 1. c) With imper- 
fections £j = 10~ 2 , £b = averaged over N av = 10 2 realiza- 
tions. 



which for j = 1, • • • , N represents a state of the chain 
where the jth spin is prepared in |1) and the other N — 1 
ones in |0). The global state of the chain at time t is 

N 

|*(t)) = cos tf|0) +sin 0e** |j) , (4) 

where |0) is the chain state with all the spins in |0) and 
where (we set H = 1) 
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The accuracy of the state transfer is determined through 
the analysis of the fidelity 

T{t, i?, <p) = (*„(0, <p)\p N (t) |*o(^, <p)) 

where pN(t) is the reduced density matrix of the iVth 
spin at time t. We consider the fidelity averaged over the 
initial state \^o) distributed uniformly over the Bloch 
sphere 3, 
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We ignore the phase of Jn as it can be gauged away by 
a proper choice of the external field. Eventually, we will 
average over different disorder realizations, that is, 



m = (Ht)h 



(7) 



|j) = |o,o, • • -0,1,0,- ■ -o) , 



(3) 



where (.)x> stands for the average over different imper- 
fection configurations. In Ref. [9j it has been shown that 



3 



after a time t n = (2n + 1)tt/4J (n integer) the state 
of the left-most spin is transferred exactly to the right- 
most spin. This is due to the fact that the Hamiltonian 
(J5J) can be viewed as that of a pseudospin S — (N — 1) /2 
that precesses in a constant magnetic field. The state 
transmission is equivalent to a t rotation of the spin. 
In order to analyze the robustness of this protocol to 




FIG. 2: Averaged fidelity at time ti as a function 
of the disorder e.j for different spin-chain lengths and 
e b = 0, iV a „ = 10 3 . From right to left N = 
10, 20, 50, 100, 150, 200, 300, 400, 500. Jnset: e} as a func- 
tion of N obtained from the condition Fiti) = 0.9 (circles) 
and J-(ti) — 0.7 (squares). Straight lines are proportional to 
jy-o.5 Here and in the following figures the logarithms are 
decimal. 



static imperfections, we model their effects by adding to 
the Hamiltonian a random perturbation both in the ex- 
change couplings and in the local variations of the mag- 
netic field. The coefficients in Eq.QJ are replaced with 
the new values 

B k ^b k , J k ^J k (l + 6 k ) 

where 5k and b k are random variables with uniform distri- 
bution in the intervals 5k £ [— £j, £j] and b k € [— £b,£b]- 
The results presented in this paper are obtained by av- 
eraging over N av different disorder realizations. 



III. STABILITY OF THE COMMUNICATION 
IN A DISORDERED CHAIN 

We numerically solve the Schrodinger equation for the 
dynamical evolution and compute the fidelity of the right- 
most spin with respect to the input state. In Fig. ^ wc 
plot typical results of this evolution both for the ideal 
case (Fig.^l) and in presence of imperfections (Fig.[Tj3.). 
Figure [T]p is the result of an average over different dis- 
order realizations. In the presence of disorder the simple 
periodicity of the fidelity oscillation is lost. Moreover, 



the maximal value of the fidelity is less than unity (it is 
reached at slightly different time intervals as compared 
to the ideal case). Thus the optimal time for state trans- 
fer should be inferred for each experimental sample. The 
original (in the ideal case) periodicity of the signal is 
recovered averaging over different disorder realizations, 
however, the maxima are progressively suppressed on in- 
creasing time. Therefore the optimal state transfer, in 
presence of imperfections, is obtained in correspondence 
of the first peak at time t\ = 7r/4J. 

In this section we concentrate on the dependence of 
the optimal fidelity [T at time t\) as a function of static 
imperfection strength and of the chain length. 
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FIG. 3: Averaged fidelity at time ti as a function of 
magnetic field disorder eb for different spin-chain lengths 
and ej = 0, N av = 10 3 . From left to right N = 
10, 20, 50, 100, 150, 200, 300, 400, 500. Jnset: e b as a func- 
tion of N obtained from the condition T{t\) = 0.9 (circles) 
and J-(ti) = 0.95 (squares). Straight lines are proportional 
to iV ' 43 . 



In Fig. |2 we report the fidelity as a function of e j for 
different chain lengths assuming, for the moment, that 
there is no disorder in the local field (eb = 0). The 
opposite situation, with disordered local magnetic field 
(eb 7^ 0) and ideal nearest neighbor interaction (ej = 0) 
is shown in Fig. |SJ These sources of disorder lead to a 
striking different behavior. While in the first case the 
error introduced by the imperfections increases with N, 
the effect of the disorder on local magnetic field decreases, 
becoming less effective on increasing the chain length. 
For completeness we show the case where both eb and 
Ej are different from zero in Fig. 21 The fact that the 
two effects are almost independent can be traced back to 
the fact that the we are working in the sector with one 
spin up. 

The behavior of the fidelity obeys a simple scaling law. 
We verified numerically that the fidelity scales as 
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T(t l ) = -{l + e 



-KjNej-K B e B /N 
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where kj ~ 0.2 and kb ~ 0.7. The constants kj (kb) 
have been obtained from the dependence, as a function 
of N, of the value £j (e B ) at which the fidelity reaches a 
given threshold value (see the insets of Figs. |21 and |3J). 

The scaling given in Eq. JSJl can be justified in the limit 
of very small disorder by means of perturbation theory. 
In the limit sjt, est <C 1, the fidelity reads 

2 N 
k=l 

2 N 

- £ -fY,( 2 ^ F »Mt)}~El(t))/3 (9) 

fe=i 

The coefficients Ck, Dk.k^Ek, and Fk,k as well as the de- 
tails of the calculation are given in the appendix The 
disorder in the local magnetic field averages out in the 
limit of infinite spin chains. In view of the little effect of 
random fields on the quantum communication over long 
chains, from now on we will consider only the effect of 
disordered exchange coupling between spins. 




FIG. 4: Average fidelity at time t\ as a function of the 
amplitudes of disorder ej, eb for a N = 50 spin network, 
No* = 100. 



The presence of spatial correlation in the disorder 
is a concrete possibility in experimental realizations of 
this protocol, as, for example, with Josephson-junction 
chains 0, 0- We model correlated disorder as follows: 
The sign of any single 5k, the error on the fcth coupling, 
is correlated with the previous one following the rule: 

SiSi—i > with probability "P, 

5i5i-i < otherwise. (10) 

The correlations introduced in Ea. (|10fl result in a per- 
fect correlation (anticorrelation) in the signs between the 
fluctuations among nearest neighbors if V = 1 (V = 0). 
Uncorrelated disorder is recovered for V = 0.5. In Fig. [S] 
the fidelity as a function of e,j (eb — 0) for different 



values of V is plotted. Notice that the fidelity decay is a 
monotonic function of V . Anticorrelated disorder is more 
dangerous than correlated one. 




FIG. 5: Fidelity at time t\ as a function of ej with N = 100, 
N av = 200 and from left to right V = 0.1, 0.25, 0.5, 0.75, 0.9. 



Within the model of disorder studied in this work, 
strong fluctuations of the exchange couplings lead to a 
degradation of the signal while the same protocol is not 
very sensitive (especially for long chains) to fluctuations 
in the local magnetic fields. 

IV. LEVEL SPACING STATISTICS 

The behavior of the fidelity is essentially dictated by 
the time dependence of the amplitude fjq defined in 
Eq. JSJ. A deeper insight of its characteristics in disor- 
dered chains can be understood by analyzing the statis- 
tics of the level spacing of the spin-chain Hamiltonian in 
presence of disorder. The level spacing statistics P(s) is 
widely used to study complex many-body systems [jjj 
and quan tum systems with classically chaotic counter- 
parts [23] in the framework of random matrix theory |16| . 
The distribution P(s)ds gives the probability that the 
energy difference between two adjacent levels (normal- 
ized to the average level spacing) belongs to the interval 
[s, s + ds]. The Hamiltonian Q is a tridiagonal matrix 
and thus it is not a random matrix, however, we will show 
that this analysis helps in understanding the behavior of 
the disordered chain. The level spacing statistics can still 
be used to characterize the crossover that static imper- 
fections induce in the spectrum of the Hamiltonian. 

The Hamiltonian H is a tridiagonal matrix with 
zero entries on the diagonal, and H^.k+i = Hk+i,k = 
\yJk{N — k) where N is the chain length and A a con- 
stant. Without any perturbation [ej = eb = 0) the 
energy levels are then equally spaced, while in presence 
of strong random perturbations (|ej| ~ 1) its eigenvalues 
are completely uncorrelated. This crossover is detected 



FIG. 6: Level spacing statistics P(s) for N = 100, es = 0, 
N av — 10 3 and different values of ej: ej = 10 -3 (circles), ej = 
2 x 10~ 2 (squares), ej = 5 x 10~ 2 (diamonds), ej = 2 X 10 _1 
(triangles up), ej = 1 (triangles down). Inset: magnification 
of the same figure around s = 1. The dashed line corresponds 
to the Poissonian Pp(s). 



by the level spacing statistics. It changes from a delta 
function to a Poisson distribution given by the formula 

Pd{s) — S(s — 1) no disorder (11) 
Pp(s) = exp(— s) strong disorder. (12) 

Figure HJ1 shows this crossover: P(s) changes from one 
limiting case to the other as a function of static imper- 
fection strength. This crossover can be quantitatively 
characterized by the parameter: 



? _ ^ \P(s) - P P (s)\d S 
V fi\P D (s)-P P (s)\ds' 



(13) 



which varies from r\ = 1 in the case of a delta function to 
r] = for a Poisson distribution 0] ■ bi Fig. we show 
the dependence of rj on the strength of the perturbation. 
The crossover starts at ej ~ 1CP 3 — 1CP 2 depending on 
the length of the chain. In the inset of Fig.[7|we report the 
dependence of the imperfection strength rj c at which the 
parameter rj reaches a given constant value (77 = 0.5, 0.8). 
The threshold r\ c drops with the spin length as 
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Thus it follows the same law found in the previous section 
regarding the fidelity of the state transfer. 

In the next section, we show that the same crossover is 
reflected by the fidelity time series with the appearance 
of a fractal behavior. 



FIG. 7: The parameter i) as a function of the strength of 
the static imperfections e.j. Different curves correspond to 
different spin chain length: from right to left N = 50, N = 
100, N = 200, N = 500. We averaged over N a v disorder 
realizations with N av = 10 4 . Inset: scaling of the parameter 
r/ c as a function of the chain length iV obtained from the 
condition r\ — 0.5 (circles) and r\ — 0.8 (squares). Straight 
lines are proportional to iV~ ' 5 . 



imperfections emerges in the time dependence of the fi- 
delity. In this section we will not look for the optimal 
time for the state transfer but rather analyze its behav- 
ior as a function of time. It appears that the time signal 
of the fidelity has a fractal behavior. In order to measure 
the fractal dimension of the sig nal we used the modi- 
fied box counting algorithm |21|. In the standard box 
counting algorithm the fractal dimension D of the signal 
is obtained by covering the data with a grid of square 
boxes of size L? . The number M(L) of boxes needed to 
cover the curve is recorded as a function of the box size 
L. The (fractal) dimension D of the curve is then defined 
as 



D = - Urn \og L M(L). 



(14) 



One finds D = 1 for a straight line, while D = 2 for a 
periodic curve. Indeed, for times much larger than the 
period, a periodic curve covers uniformly a rectangular 
region. Any given value of D in between of these integer 
values is a signal of the fractality of the curve. The mod- 
ified algorithm of Ref. [2lJ follows the same lines but uses 
rectangular boxes of size L x Aj (Aj is the largest ex- 
cursion of the curve in the region L). Then, the number 



M(L) 



(15) 



V. FRACTAL DIMENSION OF THE FIDELITY 

An interesting consequence of the modification of the 
spectrum, and hence of the fidelity, in presence of static 



is computed (the time boxes L are expressed in units of 
the exchange coupling J) . For any curve a region of box 
lengths L min < L < L max exists where M oc L D . Out- 
side this region one either finds D = 1 or D = 2: The 
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first equality [D = 1) holds for L < L m i n and it is due 
to the coarse grain artificially introduced by any numer- 
ical simulations. The second one (D = 2) is obtained for 
L > L max and it is due to the finite length of the ana- 
lyzed time series. The boundaries L m i n , L max have to be 
chosen properly for any time series. 
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log 10 (JL) 

FIG. 8: M as a function of the interval JL. A numerical fit 
gives a fractal dimension D — 1.52. Inset: Temporal evolution 
of the fidelity up to time T = 10 4 / J in the presence of disorder 
for ej = 0.26, e B = 0, N av = 1, N = 500. 

We apply the modified algorithm to the signal of the fi- 
delity for a single realization of disorder after a transient 
regime needed to reach the average value of T = 0.5: 
The inset of Fig. [S] shows the typical fluctuating signal 
we analyzed while Fig. [S] shows the numerically com- 
puted function M(V) which gives a fractal dimension 
D = 1.52. It is natural to investigate the dependence 
of the fractal dimension with the static imperfection 
strengths: the results of numerical simulations are given 
in Fig. [5] The curve changes gradually its dimension from 
fl«2 (periodic curve) to D = 1 for very large imperfec- 
tion strengths. This last result is due to the fact that for 
very large disorder the fidelity drops almost immediately 
to 0.5 corresponding to a complete loss of the initial state 
information: The fidelity remains then constant, charac- 
terized by dimension D — 1. However, the most general 
situation in presence of static imperfections is a fidelity 
with fractal dimension: defining, as before, a threshold 
of disorder strength D c at which the fidelity has a given 
fractal dimension (between two and one), we find that 
this threshold drops as 

D c ~ TV" ' 5 . 

This behavior is shown in Fig. 1101 and follows exactly the 
same scaling as the parameters rj c and e°j. 

VI. CONCLUSIONS 

We have shown that static imperfections in a mod- 
ulated spin chain destroy, above a given threshold, the 




FIG. 9: Fractal dimension D of the signal J-(t) as a function 
of the perturbation strength sj for eb = 0, N a v = 1 and, from 
right to left, TV = 100, TV = 200, TV = 500, TV = 1000. The 
error on the fractal dimension is of the order of three percent. 
For £j < 5 x 10 -2 the error on the fractal dimensions D 
increases significantly as L m i n < Lmax- 



transmission of quantum states if performed following 
the protocols presented in Ref. @- We characterize the 
effects of static imperfections by means of the fidelity 
of the state transmission. This transition is reflected in 
the changing of the level spacing statistics (from delta- 
correlated to completed uncorrelated) and in the behav- 
ior of the fidelity time series: the perfect state trans- 
fer is characterized by a periodic fidelity with integer 
fractal dimension while beyond the critical threshold it 
is described by a fractal dimension. We characterize 
these crossovers by analyzing e L j,r] c ,D c : the imperfec- 
tion strength needed to reach this value defines a criti- 
cal threshold. The three distinct critical thresholds fol- 
low the same scaling as a function of the chain length 
and imperfection strength, independently from the crit- 
ical value chosen. This common behavior reflects the 
profound changes in the quantum system induced by the 
presence of static imperfections. The threshold drops as 
the square root of the chain length: this is a behavior 
similar to the one found in Ref. [lj] in a different system 
where it was a consequence of the two body nature of the 
interactions. Here, the dependence is mainly due to the 
fact that the system is confined to the subspace of one 
excitation. The conclusion of this analysis is that it is 
possible, at least in principle, to tolerate or correct the 
errors introduced by static imperfection. 
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APPENDIX A: PERTURBATION THEORY 

We are interested in evaluating the fidelity J7J averaged 
over different disorder realizations. Equation © shows 
that it depends on the matrix element 



f N (t) = (N\e-< H+H ^\l 
= (N\e-' lHt T 



exp 



10 



-iHt 



(Al) 
|1> 



= 1 + OiHj) + 0{H 2 j) 

where T is the time ordered product, K = 1, and Hj 
is the part of the Hamiltonian that describes the static 
imperfections bk,Sk- We first consider the case where 
5k = 0, that is, only random local magnetic fields are 
present. We develop the time ordered product up to the 
second order in Hi. The first order term reads 

O(iJj) = -i [ dt(l\e lHt H ie - lHt \l) (A2) 
Jo 

N t N 

= -i^h dt(l~2\Ul(t)\ 2 ) = -iY / beC e (t), 
i=i ^° i=i 

where Uj?(t) = {l\e-' lHt \k) j^. The second order is 
given by 

0(Hj) = - f f dtdt'(l\e lHt H I e-* H ( t - t ' , >H I e-* Ht '\l) 
Jo Jo 



= -Y.Y.btbm / dtdt' l-2\ul{t)\ 
i=i m =i Jo Jo L 

N 

- 2 \u}(t')\ 2 + 4U%(t)U}(f) £ U^Wktf) 

N N 

= ~2J X! bib m D itm (t) 
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(A3) 



= 1 m=l 



The fidelity JJJ) is given by the average over different dis- 
order realization of the coefficient /./v(i) and of its mod- 
ulus square: 



.\f N (t)\ , \f N (t)\< 
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(A4) 



2 N 

l-|^(2fc[%i-C t 2 (f))/3. 

k=l 



The case for 8 k is obtained following the same steps 
and with a final result of 



2 N 
k=l 

2 N 

- e -J-J2(^4FkAt)}-E 2 k (t))/3, (A5) 



k=l 



where the coefficients E)~, F m ^ are given by 



dtRe [U}(t)U&M, 

N 

dtdt'j2iu?:(t)u™+ i (t) 



(A6) 



o Jo 

i* 



k=l 



ull^t'Wtit 1 )). 



The effects of the two different kind of perturbations in 
Eq. (|A5|) . the local magnetic bi fields and the couplings 
<5fc, are decoupled because they fluctuate independently 
from each other. 
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